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ABSTRACT 


Independent Earth orientation measurements taken by the space-geodetic techniques of lunar and 
satellite laser ranging, very long baseline interferometry, and the Global Positioning System have 
been combined using a Kalman filter. The resulting combined Earth orientation series, SPACE2016, 
consists of values and uncertainties for Universal Time, polar motion, and their rates that span from 
September 28, 1976, to June 30, 2017, at daily intervals and is available in versions with epochs 
given at either midnight or noon. The space-geodetic measurements used to generate SPACE2016 
have then been combined with optical astrometric measurements to form two additional combined 
Earth orientation series: (1) COMB2016, consisting of values and uncertainties for Universal Time, 
polar motion, and their rates that span from January 20, 1962, to June 30, 2017, at daily intervals and 
which are also available in versions with epochs given at either midnight or noon; and (2) POLE2016, 
consisting of values and uncertainties for polar motion and its rate that span from January 20, 1900, 
to June 22, 2017, at 30.4375-day intervals. 
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INTRODUCTION 


Reference series of Earth orientation parameters (EOPs) obtained by combining independent mea- 
surements of the Earth’s orientation are generated annually at the Jet Propulsion Laboratory (JPL) 
in support of tracking and navigation of interplanetary spacecraft. This report describes the gen- 
eration of the most recent such combined Earth orientation series: SPACE2016, COMB2016, and 
POLE2016. Since the procedures used to generate these most recent series are similar to those used 
to generate previous such combinations, only a brief description of their generation is given here. 
Further details regarding the approach used at JPL to annually combine independent measurements 
of the Earth’s orientation can be found in Gross (1996, 2000) and Gross et al. (1998) 


SPACE2016 


Data Sets Combined to Form SPACE2016 


SPACE2016 is a combination of independent space-geodetic measurements of the Earth’s orientation. 
Table 1 lists the space-geodetic series used in generating SPACE2016, giving their identifiers, the 
number of measurements from each series that were actually incorporated into SPACE2016, and 
the time interval spanned by those measurements. Note that the University of Texas Center for 
Space Research (UTCSR) satellite laser ranging (SLR) Universal Time (UT) values were not used 
in generating SPACE2016 due to problems associated with separating this component of the Earth’s 
orientation from the effects of unmodeled forces acting on the satellite that cause the node of its orbit 
to drift (see Gross et al. 1998, p. 217 for further discussion about this point). For similar reasons, the 
International Laser Ranging Service (ILRS) satellite laser ranging length-of-day (LOD) values have 
also not been used in generating SPACE2016. 

Since it was desirable to combine only independent measurements of the Earth’s orientation, only 
those series listed in Table 1 were used, even though other space-geodetic series are available from 
other analysis centers. When more than one series determined by the same measurement technique 
was used, care was taken to make sure that the measurements themselves were not included more 
than once. In particular, polar motion measurements from the JPL Global Positioning System (GPS) 
series were only used until the start of the International Global Navigation Satellite Systems (GNSS) 
Service (IGS) combined series EOP(IGS) 95 P 01 on January 1, 1995; polar motion measurements 
from the IGS combined EOP(IGS) 95 P 01 were then used until the start of the IGS combined 
series EOP(IGS) 95 P 02 on June 30, 1996; polar motion measurements from the IGS combined 
series EOP(IGS) 95 P 02 were then used until the start of the accumulated IGS Solution Independent 
Exchange (SINEX) combined series EOP(IGS) 00 P 03 on February 26, 2000; and polar motion 
measurements from the accumulated IGS SINEX combined series EOP(IGS) 00 P 03 were used 
thereafter. Similarly, LOD measurements from the IGS combined series EOP(IGS) 95 P 02 were 
used until it ended on February 25, 2000 with LOD measurements from the accumulated IGS SINEX 
combined series EOP(IGS) 00 P 03 used thereafter. In addition, UT measurements from the National 
Oceanic and Atmospheric Administration (NOAA) International Radio Interferometric Surveying 
(IRIS) Intensive UT1 series were used until it ended on December 31, 1994; measurements from 
the United States Naval Observatory (USNO) National Earth Orientation Service (NEOS) Intensive 
UT series were then used until it ended on December 4, 2000; and measurements from the Goddard 
Spaceflight Center (GSFC) NEOS Intensive UT1 series were used thereafter. Finally, polar motion 
measurements from the UTCSR SLR series EOP(CSR) 96 L 01 were used until it ended on February 
3, 1996, with measurements from the ILRS combined SLR series being used thereafter. 


Data Preprocessing and Treatment of Tide-Induced Rotational Variations 


The Earth orientation series listed in Table 1 were first preprocessed by removing leap seconds from 
the Universal Time (UT1) values and, when necessary, by correcting the UT1 values to be consistent 
with the extended definition of Greenwich Sidereal Time (GST) as adopted by the International Earth 
Rotation and Reference Systems Service (ERS; IERS 1997, p. 149). Since most of the series listed 
in Table 1 were already consistent with the extended definition of GST, this correction needed to be 
applied to only the NOAA IRIS Intensive UT1 series. 


Table 1. Data Sets Combined to Form SPACE2016* 


Data Analysis 


Data Set Name Type Center Reference Data Span Number 
LLR (JPL17M01; VOL, UTO) 

McDonald Cluster LLR JPL Williams et al. (2012) Oct. 5, 1976, to Jan. 5, 2012 742 

OCA LLR- JPL Williams et al. (2012) Apr. 7, 1984, to Jun. 25, 2017 1485 

Haleakala LLR- JPL Williams et al. (2012) Feb. 10, 1985, to Aug. 11, 1990 70 

Apache Point LLR- JPL Williams et al. (2012) Jun. 4, 2006, to Oct. 15, 2016 519 
UTCSR (CSR96L01) 

LAGEOS Polar Motion SLR UTCSR Eanes & Watkins (1996) Sep. 27, 1976, to Feb. 4, 1996 2218 
ILRSA (11Aug17; Polar motion) 

ILRS Primary Combination SLR ASI-CGS Sciaretta et al. (2010) Feb. 5, 1996, to Jun. 29, 2017 7166 
DSN (JPL97RO1; T, V) 

California-Spain Cluster VLBI JPL Steppe et al. (1997) Nov. 26, 1979, to Jun. 3, 2017 1222 

California-Australia Cluster VLBI JPL Steppe et al. (1997) Oct. 28, 1978, to Sep. 30, 1997 698 
NOAA (NOAA95R02) 

IRIS Intensive UT1 VLBI NOAA Ray et al. (1995) Apr. 2, 1984, to Dec. 31, 1994 2393 
USNO (N9903) 

NEOS Intensive UT1 VLBI USNO Eubanks et al. (1999) Jan. 4, 1995, to Dec. 4, 2000 1497 
NASA/GSFC (GSF2016a) 

NEOS Intensive UT1 VLBI GSFC NASA/GSFC, VLBI Group (2016) Dec. 6, 2000, to Jun. 29, 2017 5375 
NASA/GSFC (GSF2016a) 

Multibaseline VLBI GSFC NASA/GSFC, VLBI Group (2016) Apr. 12, 1980, to Jun. 30, 2017 5344 
NASA/GSFC (GSFC1122) 

Westfort-Fort Davis VLBI GSFC Gordon et al. (1999) Jun. 25, 1981, to Jan. 1, 1984 105 

Westfort-Mojave VLBI GSFC Gordon et al. (1999) Mar. 21, 1985, to Aug. 6, 1990 18 


GPS (21APR04; Polar motion) 
Post-processed Flinn Analysis GPS JPL Heflin et al. (2004) Jun. 10, 1992, to Dec. 31, 1994 817 


GPS (IGS95P01; Polar motion) 
IGS Final Combined (ITRF93) GPS NRCan Kouba & Mireault (1997) Jan. 1, 1995, to Jun. 29, 1996 546 


GPS (IGS95P02; Polar motion) 
IGS Final Combined (IGS08) GPS CODE Mireault & Kouba (2000) Jun. 30, 1996, to Feb. 26, 2000 1337 


GPS (IGSOOP03; Polar motion) 
IGS SINEX Combined (IGSO8) GPS IGN Rebischung & Garayt (2013) Feb. 27, 2000, to Jun. 29, 2017 6321 


GPS (IGS95P02; LOD) 
IGS Final Combined (IGS08) GPS CODE Mireault & Kouba (2000) Feb. 23, 1997, to Feb. 26, 2000 1098 


GPS (IGSOOP03; LOD) 
IGS SINEX Combined (IGS08) GPS IGN Rebischung & Garayt (2013) — Feb. 27, 2000, to Jun. 29, 2017 6330 


* LLR, lunar laser ranging; JPL, Jet Propulsion Laboratory; VOL, variation of latitude; UT, Universal Time; OCA, Observatoire de la 
Céte d’Azur; UTCSR, University of Texas Center for Space Research; LAGEOS, Laser Geodynamic Satellite; SLR, satellite laser rang- 
ing; ILRS, International Laser Ranging Service; ASI, Agenzia Spaziale Italiana; CGS, Centro di Geodesia Spaziale; DSN, Deep Space 
Network; T, transverse; V, vertical; VLBI, very long baseline interferometry; NOAA, National Oceanic and Atmospheric Administra- 
tion; IRIS, International Radio Interferometric Surveying; USNO, United States Naval Observatory; NEOS, National Earth Orientation 
Service; NASA, National Aeronautics and Space Administration; GSFC, Goddard Space Flight Center; GPS, Global Positioning Sys- 
tem; IGS, International Global Navigation Satellite System (GNSS) Service; ITRF, International Terrestrial Reference Frame; NRCan, 
Natural Resources Canada; CODE, Center for Orbit Determination in Europe; IGN, Institut National de l’ Information Géographique et 
Forestiére; SINEX, solution independent exchange; LOD, length of day. 
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Table 2. Changes to the IGS Reference Frame. 


Date of Change Previous Frame New Frame _ Reference 


June 30, 1996 ITRF93 ITRF94 IGS Analysis Center Coordinator (2012) 
March 1, 1998 ITRF94 ITRF96 Kouba et al. (1998) 

August 1, 1999 ITRF96 ITRF97 Springer (1999) 

February 27,2000 ITRF97 IGS97 Ferland (2000) 

December 2, 2001 IGS97 IGSO00O Weber (2001) 

January 11, 2004 IGSO00 IGb00 Ferland (2003) 

November 5, 2006 IGb00 IGSO5 Ferland (2006) 

April 17, 2011 IGSO5 IGS08 Rebischung et al. (2011, 2012) 

October 7, 2012 IGS08 IGb08 Rebischung (2012) 

January 29, 2017 IGbO08 IGS14 Rebischung (2016) 


Rotational variations caused by solid-Earth and ocean tides were also removed from the UT1 
values. The effect of the solid-Earth tides was removed by using the model of Yoder et al. (1981); the 
model of Kantha et al. (1998) was used to remove the effect upon UT1 of the ocean tides at the Mf, 
Mf’, and Mm tidal frequencies. Since the Yoder et al. (1981) model already includes a contribution 
from the equilibrium ocean tides, only the Kantha et al. (1998) oceanic corrections to the Yoder 
et al. (1981) model were actually removed. Also note that the Kantha et al. (1998) model was used 
to remove the effect of ocean tides on only UTI, not on polar motion. Ocean-tide-induced polar 
motion variations were not removed from any of the polar motion observations. Finally, the only 
Earth orientation series listed in Table 1 that includes the effects of semidiurnal and diurnal ocean 
tides on the Earth’s orientation is the NOAA IRIS Intensive UTI series. This series included these 
effects by adding to the released UT1 values the model of Herring (1993, also see Herring & Dong 
1994). Hence, the same Herring (1993) model was used to remove them. 

From time to time, the IGS has changed the terrestrial reference frame used to express its prod- 
ucts, usually related to updates of the International Terrestrial Reference Frame (ITRF). A chronology 
of the relevant IGS changes is given in Table 2. Such changes in reference frame potentially introduce 
discontinuities into the IGS combined Earth orientation series. However, the IGS Final combined se- 
ries EOP(IGS) 95 P 02 and the accumulated IGS SINEX combined series EOPUGS) 00 P 03 used 
here have had each of these segments aligned to the same IGS reference frame. Thus, to within the 
uncertainty in determining the corrections required to align each segment, the IGS combined series 
EOP(IGS) 95 P 02 and EOP(IGS) 00 P 03 used here should be reasonably consistent with each other. 
They were, therefore, concatenated with one common set of bias-rate corrections being determined 
for them, as described below. 


Adjustments Made to Space-Geodetic Series Prior to Combination 


Prior to combining the series listed in Table 1 to form SPACE2016, series-specific corrections were 
applied for bias and rate, the stated uncertainties of the measurements were adjusted by multiplying 
them by series-specific scale factors, and outlying data points were deleted. Values for the bias-rate 
corrections and uncertainty scale factors were determined by an iterative procedure wherein each 
series was compared, in turn, to a combination of all others. In order to minimize interpolation error 
(see Gross et al., 1998, pp. 223-225), the comparison of each series to its reference combination was 
done at the epochs of the measurements of that series by generating its reference combination using a 


Kalman filter that interpolates to (and prints its EOP estimates at) the exact epochs of those measure- 
ments. Also, both the bias-rate corrections and the uncertainty scale factors for all components of 
a given series were determined simultaneously in a multivariate approach using nonlinear weighted 
least squares. Using a multivariate approach allows the correlations between the components to 
be taken into account when determining the bias-rate corrections and uncertainty scale factors (see 
Gross et al., 1998, pp. 22). 

All the series listed in Table 1 except for the Westford—Mojave and Westford—Fort Davis base- 
lines of GSFC VLBI were included in the iterative procedure. Details of the iterative procedure are 
described in Gross (1996, 2000) and Gross et al. (1998) and will not be repeated here. These details 
include: 


1. Use of a reference series, SPACE2015 (Ratcliff & Gross, 2016), for initial bias-rate alignment; 
2. Analysis of each data type in its natural reference frame; 


3. Clustering of the McDonald LLR stations and, separately, the DSN very long baseline interferometry 
(VLBI) stations in California, Australia, and Spain; 


4. Initial convergence on values for the series-specific bias-rate corrections and uncertainty scale factors 
prior to the start of outlier detection and deletion; and 


5. Final convergence on these values after detecting and deleting all data outliers. 


At the end of the iterative procedure, relative bias-rate corrections have been determined that 
make the series agree with each other in bias and rate; uncertainty scale factors have been determined 
that make the residual of each series have a reduced chi-square near one when differenced with a 
combination of all others; and outlying data points (those with residual values greater than four times 
their adjusted uncertainties) have been deleted. A total of 121 data points, or about 0.26% of those 
combined, were thus deleted. 

For GSFC NEOS Intensive and GSFC Multibaseline VLBI, the initial alignment to SPACE2015 
was performed individually for each baseline containing at least 90 data points and spanning more 
than | year. Baselines containing less than 1 year of data or fewer than 90 data points were aligned 
as a group to SPACE2015. Subsequent steps of the iterative procedure were then conducted indepen- 
dently on the individual baseline groups. 

Starting April 17, 2011 the IGS changed their procedure for estimating uncertainty in GPS EOP 
values. This introduced a change in magnitude of the polar motion X (PMX), polar motion Y (PMY), 
and LOD uncertainties of approximately a factor of 6. To account for this sudden change, uncertainty 
scale factor values after April 17, 2011 were estimated separately from those of earlier epochs using a 
trial reference series. These trial scale factors were applied to the original IGS GPS data, which were 
then used in the iterative procedure described above to form a final reference series. The resulting 
scale factors shown in Table 3 for the IGSOOPO3 data set include the contribution of both the trial 
scale factors and the overall iteratively determined scale factor. 

Since the GSFC VLBI data from the Westford—Mojave and Westford—Fort Davis baselines con- 
tain relatively few measurements over a short time span, they were not included in the iterative pro- 
cedure. Instead, relative bias-rate corrections for these VLBI series were determined that make them 
agree in bias and rate with an intermediate combination that was obtained by combining all the other 
series after applying to them the above iteratively determined bias-rate corrections and uncertainty 
scale factors. Uncertainty scale factors were determined that make the VLBI residuals have a reduced 
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Table 3. Adjustments to Space-Geodetic Series* 


Data Set Name Bias, mas Rate, mas/yr Uncertainty Scale Factor 
LLR (JPL17M01) VOL UTO VOL UTO VOL UTO 
McDonald Cluster — 0.1366 + 0.0875 — 0.0725 + 0.0880 0.0558 + 0.0204 0.0211 + 0.0181 5031 1.3584 
OCA 0.0002 + 0.0287 — 0.0683 + 0.0254 — 0.0008 + 0.0022 0.0041 + 0.0018 2.1017 1.7338 
Haleakala — 1.6400 + 1.1572 — 1.1084 + 0.7210 — 0.2060 + 0.2252. — 0.1551 40.1512 2.7426 2.4148 
Apache Point — 0.0568 + 0.1764 — 0.3328 + 0.1808 0.0018 + 0.0097 0.0205 + 0.0096 9362 1.9469 
DSN (JPL97RO1) T Vv T Vv T Vv 
California—Spain Cluster — 0.7014 40.0220 — 0.3096 + 0.0553 0.1777 + 0.0012 0.1810 + 0.0034 2.0742 1.4524 
California—Australia Cluster 0.5747 + 0.0184 0.4181 + 0.0526 0.0789 + 0.0076 0.0422 + 0.0236 4324 1.1667 
NASA GSFC (1122) Ty Vv T. Vv T Vv 
Westford—Fort Davis 9.6393 + 3.8370 4.1995 + 5.9178 0.9082 + 0.3858 0.4439 + 0.5925 .2651 0.9102 
Westford—Mojave 0.6529 + 1.0902 — 0.3922 + 1.6838 0.0008 0.0020 2.4784 0.7250 
GSFC (GSF2016a)° PMX PMY UTI PMX PMY UTI PMX PMY UTI 
G3GcKuRi 0.2519 + 0.1031 —0.1353+0.1051 0.1733 40.1686 0.1445 + 0.0625 — 0.0635 + 0.0652 — 0.1481 40.1023 1.3819 1.4012 1.9307 
HrRiwfWz 0.5458 + 0.1717 — 0.6679 0.1780 0.3190 + 0.2828 + ~—- 0.0699 + 0.0292 — 0.0712 + 0.0299 =: 0.0070 + 0.0481 ~=—s:1.1825. 1.2386 ~—:1.7094 
MoRiWfWz 0.1740 + 0.1186 — 0.5748 +0.1146 0.0639 40.1442 — 0.0261 + 0.0506 — 0.1457 + 0.0482 — 0.0341+40.0656 1.4257 1.3425 0.5677 
Other Networks 0.0169 + 0.0051 —0.1657+0.0046 0.2606 + 0.0088 —0.0011+0.0003 0.0026 40.0003 —0.0082+0.0005 1.6748 1.6934 1.7588 
NOAA (95R02) PMX PMY UTI PMX PMY UTI PMX PMY UTI 
IRIS Intensive ee 0.225540.0226 0 © wae 0.0544 40.0068 = ----- === 0175 
USNO (N9903) PMX PMY UTI PMX PMY UTI PMX PMY UTI 
NEOS Intensive ee — 0.0228 40.0396 wee 0.0278 + 0.0073 ----- = 3619 
GSFC (GSF2016a)° PMX PMY UTI PMX PMY UTI PMX PMY UTI 
KKNy ree eee 0.060640.494. 000 wre wee 0.001140.0257 = ----- — ----- .2650 
KkWz 0.192440.0183 9:2 ae —0.005540.0011 9 ---- = 5201 
TsWz —0.062240.0442 00 we 0.0054 40.0025.  ----- == 7362 
NyTsWz —0.102840.1219 wre eee 0.0071 40.0065 9 ----- —— ----- 5186 
KkSvWz 0.246040.1012 0 ©2220 ee — 0.0074 40.0054 = ----- = .7033 
Other Networks eee 0.3540£0.1350 4:42 wee —0.016740.0061 9 ----- == 9186 
UTCSR (96L01) PMX PMY UTI PMX PMY UTI PMX PMY UTI 
LAGEOS — 0.305440.0101  0.6723+0.0090  —_----- 0.0484 + 0.0036 0.15304 0.0033 ss ----- 0.7329 0.7166 ----- 
ILRSA (11AUG17) PMX PMY LOD PMX PMY LOD PMX PMY LOD 
Primary Comb. 0.0712 + 0.0050 —0.0675+0.0045  —----- — 0.0066 + 0.0003 0.0003 + 0.0003. —si----- 3.3332 3.1106 9 ----- 
GPS (21APR04) PMX PMY LOD PMX PMY LOD PMX PMY LOD 
JPL Post-Flinn — 0.22004 0.0157 0.09394 0.0130  — ----- — 0.0091 + 0.0145 —0.0268+0.0119 —----- 1.3050 0.9243 ----- 
GPS (IGS95P01) PMX PMY LOD PMX PMY LOD PMX PMY LOD 
Final Combined — 0.19244 0.0548 0.206640.0541 — ----- 0.2333 40.0191  0.303240.0190  — ----- 1.2145 0.9882 _— ----- 
GPS (IGS95P02) PMX PMY LOD PMX PMY LOD PMX PMY LOD 
Final Combined 0.0565 + 0.0039 — 0.0250 + 0.0035 — 0.0067 + 0.0088 — 0.0031 +0.0002  0.0022+0.0002 0.00074 0.0005 2.0863 1.6430 1.8267 
GPS (IGSOOP03) PMX PMY LOD PMX PMY LOD PMX PMY LOD 
before 4/17/2011 0.0565 + 0.0039 — 0.0250 + 0.0035 — 0.0067 + 0.0088 — 0.0031 +0.0002  0.0022+0.0002 0.0007+0.0005 2.0863 1.6430 1.8267 
after 4/17/2011 0.0565 + 0.0039 — 0.0250 + 0.0035 — 0.0067 + 0.0088 — 0.0031+0.0002  0.0022+0.0002  0.0007+0.0005 12.9274 8.9646 11.8671 


a) Reference date for bias-rate adjustment is 1993.0. See Table 1 footnotes. mas, milliarcseconds; PMX, polar motion X; PMY, polar motion Y. 
b) GSFC VLBI Station IDs from ftp: //cddis.gsfc.nasa.gov/vlbi/ivscontrol/ns-—codes.txt 


G3 = NRAO 85-3; Ge = Gilmore Creek; Hr = HRAS at Fort Davis; Kk = Kokee Park; Ku = 9m at Kokee Park; Mo = 12m at Goldstone; Ny = Ny Alesund, 


Svalbard, Norway; Ri = Richmond, FL; Sv = Svetloe, Russia; Ts = Tsukuba, Japan; Wf = Westford, MA; Wz = Wettzell, Germany. 


chi-square near one when differenced with the intermediate combination. The adjusted VLBI data 
were then combined with the iteratively aligned data for subsequent alignment with ITRF20014. 

Finally, each series was placed within a particular IERS reference frame by applying to it an 
additional bias-rate correction that is common to all the series. This additional correction was de- 
termined by first combining all the series after applying to them the relative bias-rate corrections 
and uncertainty scale factors that had been determined for them as described above. This interme- 
diate combination was then compared to the Earth orientation parameters of ITRF2014 (Altamimi 
et al., 2016a,b). Only the ITRF2014 EOP values during the interval January 1, 1994 — February 14, 
2015 were used for the comparison. Blatantly large outliers in this interval were removed from the 
ITR2014 series prior to the comparison. The additional bias-rate correction determined from this 
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comparison was then applied to the intermediate combination in order to make it, and therefore each 
individual series, agree in bias and rate with the ITRF EOP series. 

The total bias-rate corrections and uncertainty scale factors determined by the procedures out- 
lined above are given in Table 3. The values of the bias-rate corrections given in Table 3 are the 
sum of: (1) the initial corrections applied to each series to align them to a common reference, (2) all 
the incremental corrections applied during the iterative procedure, and (3) the additional, common 
correction applied in order to place each series within that particular reference frame defined by the 
ITF2014 Earth orientation series (Altamimi et al., 2016b). The values of the uncertainty scale factors 
given in Table 3 are the products of all the incremental scale factors determined during the iterative 
procedure. The uncertainties of the bias-rate corrections given in Table 3 are the lo standard errors of 
the incremental bias-rate corrections determined during the last iteration. There are no bias-rate en- 
tries in Table 3 for components that were either not used (e.g., the UTCSR SLR UT1 component) or 
not available (e.g., the NOAA IRIS Intensive polar motion components). Note that the same common 
ITRF2014 rate correction is applied to all the data sets, including those such as the Westford—Mojave 
single baseline VLBI series, for which no relative rate correction could be determined. Therefore, the 
rate corrections given in Table 3 for those series for which no relative rate corrections could be deter- 
mined are simply the common ITRF2014 rate correction, but given in the natural reference frame of 
that series. In these cases, uncertainties for the rate corrections are not given. Also note that the bias- 
rate corrections in Table 3 for the IGS combined series EOPUIGS) 95 P 02 and EOP(IGS) 00 P 03 
are the same. Since these series were initially given within the same IGS reference frame, the series 
were merged, and a common bias-rate correction was determined for them. However, uncertainty 
scale factors for these series do differ after April 17, 2011, for the reasons described on page 5. 


Combined EOP Series: SPACE2016 


A Kalman filter was used to combine the series listed in Table 1 after the bias-rate corrections and 
uncertainty scale factors listed in Table 3 had been applied to them. The resulting combined Earth 
orientation series, SPACE2016, consists of values (Figure 1) and lo standard errors (Figure 2) for 
polar motion, Universal Time, and their rates spanning September 28, 1976, to June 30, 2017, at 
daily intervals; and it is available in versions for which the epochs are given at either midnight or 
noon. Leap seconds have been restored to UT1, and the effects of the solid Earth and ocean tides 
have been added back to the UT1 values using the same models for these effects that were originally 
used to remove them from the raw series: namely, the Yoder et al. (1981) model for the solid Earth 
tides and the Kantha et al. (1998) model for the ocean tides. However, semidiurnal and diurnal ocean 
tidal terms have not been added to and are therefore not included in the SPACE2016 UTI values. 

Figure 3 is a plot of the difference between the SPACE2016 polar motion, UT1, and LOD values 
and those of the ITRF2014 Earth orientation series. Tabulated root-mean-square (rms) differences 
between SPACE2016 and ITRF2014 EOPs are inset. SPACE2016 is very consistent with ITRF2014 
EOPs, especially after January 1, 2000, where the rms of their differences are 0.042 milliarcseconds 
(mas) in the x-component of polar motion, 0.037 mas in the y-component, 0.088 milliseconds (ms) 
in UTI, and 0.069 ms in LOD. Prior to 2000, the differences exhibit greater variability. The bla- 
tant outliers that were removed from the ITRF2014 series for determining the bias-rate corrections 
have been retained for this comparison and can be seen in panel (c) near the years 1996 and 2000. 
SPACE2016 is also very consistent IERS Bulletin A (Stamatakos et al., 2015) as shown in Figure 4 
where rms differences after January 1, 2000 are 0.047 mas in the x-component, 0.043 mas in the 
y-component, 0.011 ms in UT1, and 0.012 ms in LOD. 
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COMB2016 


COMB2016 extends SPACE2016 by additionally incorporating the optical astrometric polar motion 
and UT1 series that was determined at the Bureau International de |’ Heure (BIH) from an analysis of 
time and latitude observations by Li (1985, also see Li & Feissel 1986). This BIH optical astrometric 
series consists of values and uncertainties for polar motion and UT1 that span from January 5.0, 1962, 
to December 31.0, 1981, at 5-day intervals. Also incorporated in COMB2016 are the McDonald 
Observatory LLR data that span from April 15, 1970 to September 13, 1976 that were not included 
in SPACE2016. These additional LLR points were adjusted in bias, rate, and uncertainty as per the 
McDonald Cluster in Table 3. 


Data Preprocessing and Treatment of Tide-Induced Rotational Variations 


The BIH optical astrometric series was first preprocessed by removing leap seconds from the UT1 
values and by correcting the UT1 values to be consistent with the extended definition of GST, as 
adopted by the IERS (IERS, 1997, p. 149). Rotational variations caused by solid Earth and ocean 
tides were also removed from the UT1 values. The same models that were used to remove the tidal 
effects from the series combined to form SPACE2016 were also used to remove them from the BIH 
series: namely, the Yoder et al. (1981) model for the solid Earth tides and the Kantha et al. (1998) 
model for the Mf, Mf’, and Mm ocean tides. However, since the BIH UT1 measurements represent 
an average value over a 5-day-long observation window, and since 5 days is a substantial fraction of 
the monthly and shorter-period tides, the amplitudes of these solid Earth and ocean tidal terms were 
attenuated prior to their removal from the BIH UT1 measurements. (See Gross 1996, p. 8735 and 
Gross et al. 1998, pp. 226-227 for further discussion about this point.) 


Adjustments Made to BIH Series Prior to Combination 


The preprocessed BIH optical astrometric series was combined with the space-geodetic series that 
comprise SPACE2016 after first: 


1. correcting the BIH series to have the same bias, rate, annual terms, and semiannual terms as 
SPACE2016; 


2. applying a constant multiplicative scale factor to the measurement uncertainties of the BIH 
series so that its residual, when differenced with SPACE2016, had a reduced chi-square of one; 
and 


3. deleting those data points, if any, for which residual values were greater than four times their 
adjusted uncertainties. 


The above adjustments were determined in a multivariate approach as was done for the series com- 
bined to form SPACE2016. The procedure used to determine these bias-rate and seasonal term 
corrections, uncertainty scale factors, and outlying data points has been described before (Gross, 
1996, pp. 8735-8738) and will not be repeated here. The annual and semiannual terms of the BIH 
series were adjusted in order to correct for systematic, seasonally varying effects that are known to 
be present in optical astrometric measurements. Since the values of both the BIH series and the 
SPACE2016 series are given at midnight, interpolation error (see Gross et al., 1998, pp. 223-225) 
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is automatically minimized when differencing these two series for the purpose of determining the 
adjustments to be made to the BIH series. Tables 4 and 5 give the resulting uncertainty scale factors 
and values and lo standard errors of the corrections to the bias, rate, annual terms, and semiannual 
terms thus determined for the BIH series. When determining these uncertainty scale factors and the 
corrections to the bias, rate, and seasonal terms, no outlying data points were detected. 


Combined EOP Series: COMB2016 


A Kalman filter was used to combine the BIH series with the adjusted space-geodetic series that 
comprise SPACE2016 after first applying to the BIH series the uncertainty scale factors and correc- 
tions to the bias, rate, annual, and semiannual terms that are given in Tables 4 and 5. The resulting 
combined Earth orientation series, COMB2016, consists of values (Figure 5) and 1o standard errors 
(Figure 6) for polar motion, Universal Time, and their rates that span from January 20, 1962, to June 
30, 2017, at daily intervals and is available in versions with epochs given at either midnight or noon. 
Leap seconds have been restored to UT1, and the effects of the solid Earth and ocean tides have been 
added back to the UT1 values using the same models for these effects that were originally used to 
remove them: namely, the Yoder et al. (1981) model for the solid Earth tides and the Kantha et al. 
(1998) model for the ocean tides. The full amplitude (i.e., no tidal terms attenuated) of the effects 
of the solid Earth and ocean tides at the epoch of the time tag were added back to the UTI values. 
Semidiurnal and diurnal ocean tidal terms have not been added to and are, therefore, not included in 
the COMB2016 UTI values. 


Table 4. Adjustments to Bias, Rate, and Stated Uncertainty of Optical Astrometric Series* 


Data Set Bias, mas Rate, mas/yr Uncertainty Scale Factor 
PMX PMY UT1 PMX PMY UT1 PMX PMY UT1 

BIH — 1.4879 40.8525 — 1.0485 40.6985 — 8.4242+41.0757 1.0449 + 0.5258 1.1986 + 0.4276 5.6004 40.6732 1.8119 1.6758 1.9388 

ILS — 49.4984 42.3558 —3.560841.6635 — ----- 0.4154+0.4814 —0.0091+0.3392 ~~ ----- 2.1577 1.4982 — ----- 


*Reference date for bias-rate adjustment of BIH series is 1980.0. Reference date for bias-rate adjustment of ILS series is 1970.0. 


Table 5. Adjustment to Annual and Semiannual Terms of Optical Astrometric Series* 


Data Set Coefficent of Sine Term, mas Coefficient of Cosine Term, mas 
PMX PMY UTI PMX PMY UTI 
BIH 
Annual — 5.3814+ 1.0501 — 7.8387 + 0.8814 4.9808 + 1.3881 — 2.6812 + 1.1211 9.4925 +0.9222 — 0.8861 + 1.4432 
Semiannual 2.2666 + 1.0722 —0.5840+0.8964 — 0.6147 + 1.4146 0.8682 + 1.0952 1.1277 + 0.9033 1.9790 + 1.4145 
ILS 
Annual — 2.1612 + 3.2862 7.7886 42.3123  --- 8.2447 + 3.3222 — 104683423412 = --- 
Semiannual 0.2444 + 3.2947 8.708442.3194 = —— ----- 2.2006 + 3.3108 1.8224+42.3320 — ----- 


*Reference date for adjustment of annual and semiannual terms of BIH series is 1980.0. 
Reference date for adjustment of annual and semiannual terms of ILS series is 1970.0. 
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POLE2016 


No optical astrometric observations taken at the stations of the International Latitude Service (ILS) 
were used when creating the BIH optical astrometric series used in COMB2016 (Li, 1985; Li & Feis- 
sel, 1986). The ILS polar motion measurements (Yumi & Yokoyama, 1980), which are based solely 
upon latitude observations made at the ILS stations are, therefore, independent of those comprising 
COMB2016 and have therefore been combined with them to form POLE2016. Being based solely 
upon latitude observations, the ILS series contains no UT1 measurements, but consists solely of polar 
motion measurements that span 1899.8—1979.0 at monthly intervals. Although no uncertainties are 
given with the individual polar motion values, the precision with which they have been determined is 
estimated to be 10-20 mas (Yumi & Yokoyama, 1980, p. 27). An initial uncertainty of 15 mas was, 
therefore, assigned to each of the ILS polar motion values. Since this assigned measurement uncer- 
tainty will be adjusted later, its initial value is arbitrary as long as it is not zero and serves merely as 
an a priori estimate to be used in the series adjustment procedure described below. 
The ILS series was combined with COMB2016 to form POLE2016 after: 


1. correcting the ILS series to have the same bias, rate, annual terms, and semiannual terms as 
COMB2016; 


2. applying a constant multiplicative scale factor to the measurement uncertainties of the ILS 
series so that its residual, when differenced with COMB2016, had a reduced chi-square of one; 
and 


3. deleting those data points, if any, for which residual values were greater than four times their 
adjusted uncertainties. 


These adjustments were also determined using a multivariate approach for fitting a bias, a rate, 
and these seasonal terms to the difference of the ILS series with COMB2016 during 1962.0 to 1979.0. 
The measurement uncertainties of the ILS polar motion values were adjusted by determining and 
applying a scale factor that made the residual of this fit have a reduced chi-square of one. During 
this procedure to determine uncertainty scale factors and bias, rate, and seasonal term corrections, 
no outlying ILS data points were deleted since no data points had residual values greater than four 
times their adjusted uncertainties. Tables 4 and 5 (in the COMB2016 section) also give the resulting 
uncertainty scale factors and values and lo standard errors of the corrections to the bias, rate, annual 
terms, and semiannual terms thus determined for the ILS series. 

A Kalman filter was then used to combine the ILS series with the adjusted BIH and space-geodetic 
series that comprise COMB2016, after applying to the ILS series the uncertainty scale factors and 
corrections to the bias, rate, annual, and semiannual terms that are given in Tables 4 and 5. The 
resulting combined Earth orientation series, POLE2016, consists of values (Figure 7a and 7b) and 
lo standard errors (Figure 7c and 7d) for polar motion and its rate that span from January 20, 1900, 
to June 22, 2017, at 30.4375-day intervals. 
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DISCUSSION 


The Kalman filter that was used here to combine Earth-orientation measurements contains a stochas- 
tic model of the process that is used to propagate the state vector and its associated state covari- 
ance matrix forward in time to the epoch of the next measurement. For polar motion excitation, the 
stochastic model includes a random-walk term with equal noise forcing both the x- and y-components 
of excitation and originally having a white-noise power spectral density of 246.4 mas?/day (Morabito 
et al., 1988). This level of polar motion excitation process noise in the Kalman filter was increased to 
739.2 mas”/day in order to better match the observed spectrum of polar motion excitation. Increasing 
the excitation process noise reduces the level of smoothing applied to the polar motion components of 
the propagated state vector and increases the covariance associated with those components. Thus, the 
SPACE2016, COMB2016, and POLE2016 polar motion and polar-motion excitation values are not as 
heavily smoothed as were those of SPACE2003, COMB2003, POLE2003, and earlier combinations 
produced at JPL; and the uncertainties assigned to the SPACE2016, COMB2016, and POLE2016 
polar motion and polar-motion excitation values are somewhat larger. 

The LOD components of the IGS Final combined EOP(IGS) 95 P 02 and EOP (IGS) 00 P 03 
series were incorporated into JPL’s combinations starting with SPACE2007 and COMB2007. In 
general, GPS LOD estimates are contaminated by orbital artifacts that must be removed prior to 
combination with other LOD or UT1 measurements (Chin et al., 2009). When determining the IGS 
Final combined LOD series, however, the IGS applies constraints that have the effect of remov- 
ing these orbital artifacts. Therefore the IGS Final combined LOD series was incorporated into the 
Kalman filter as a true LOD measurement for which uncertainties are assumed to be “white” (i.e., 
random, uncorrelated in time, and with a flat power spectral density). Including the daily IGS Final 
combined LOD series in JPL’s combinations helps to compensate for less frequent UT1 measure- 
ments and was found to improve the agreement of JPL’s combined LOD estimates with independent 
atmospheric and oceanic angular momentum values, especially at the highest frequencies. 

Since a Kalman filter has been used to generate SPACE2016, COMB2016, and POLE2016, the 
resulting polar motion and UT1 values are smoothed to a degree that depends upon both the spacing 
between the measurements being combined (which determines how far the state vector and state co- 
variance matrix must be propagated forward in time) and the uncertainties that have been assigned to 
the measurements. Since improvements to the observing systems, both in the hardware and software 
and in the number of systems, have led to increasingly precise determinations of the Earth’s orien- 
tation, and since the time resolution of the measurements has generally increased in concert with 
the measurement precision, the degree of smoothing applied to the SPACE2016, COMB2016, and 
POLE2016 values is a function of time, with the earlier values being more heavily smoothed than the 
more recent values. 

Daily EOP values are reported in SPACE2016 since the NOAA IRIS and GSFC NEOS Intensive 
UTI values are given at daily intervals, as are the GPS and ILRSA combined SLR values (although 
gaps exist in each of these data sets). However, prior to the start of these data sets, the measurements 
combined to form SPACE2016 are given less frequently; therefore, the Kalman filter used to combine 
these measurements also interpolates them in order to produce a series of equally spaced values. In 
order to be consistent with SPACE2016, daily EOP values are also reported in COMB2016 even 
though the BIH optical astrometric series used in COMB2016 is given at 5-day intervals. Thus, 
SPACE2016, COMB2016, and POLE2016 are equally spaced series of smoothed, interpolated Earth- 
orientation parameters. 
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The combined Earth-orientation series SPACE2016, COMB2016, and POLE2016 are available 
from JPL’s Geodynamics and Space Geodesy Group via anonymous ftp at: 
<ftp://keof.jpl.nasa.gov/keof/combinations/2016> and upon request from the 
authors: Todd.Ratcliff@jpl.nasa.gov or Richard.Gross @jpl.nasa.gov. 
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LAGEOS 
LLR 
LOD 


Acronyms and Terms 


Agenzia Spaziale Italiana 
Bureau International de 1’ Heure 


Centro di Geodesia Spaziale 

Center for Orbit Determination in Europe 

extends the SPACE series by additionally incorporating BIH optical astro- 
metric measurements of polar motion and UT1 


Deep Space Network 
Earth orientation parameter 


IVS station code for 85-3 antenna at NRAO Green Bank, WV 
IVS station code for Gilmore Creek, Fairbanks, AK 

Global Navigation Satellite System 

Global Positioning System 

Goddard Space Flight Center 

Greenwich Sidereal Time 


IVS station code for HRAS 
Harvard Radio Astronomy Station at Ft. Davis, TX 


Identification 

International Earth Rotation and Reference Systems Service 
Institut National de l’Information Géographique et Forestiére 
International Global Navigation Satellite Systems (GNSS) Service 
International Laser Ranging Service 

ILRS Primary Combination Solution 

International Latitude Service 

International Radio Interferometric Surveying 

International Terrestrial Reference Frame 

International VLBI Service for Geodesy & Astrometry 


Jet Propulsion Laboratory 


IVS station code for Kokee Park, Kauai, HI 
IVS station code for 9-m antenna at Kokee Park, Kauai, HI 


Laser Geodynamics Satellites 


lunar laser ranging 
length-of-day 
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milliarcsecond 
IVS station code for 12-m antenna at Goldstone, CA 
millisecond 


National Aeronautics and Space Administration 
National Earth Orientation Service 

National Oceanic and Atmospheric Administration 
National Radio Astronomy Observatory 

Natural Resources Canada 

IVS station code for Ny Alesund, Svalbard, Norway 


Observatoire de la Céte d’ Azur 


polar motion X 

polar motion Y 

extends the COMB series by additionally incorporating ILS optical astro- 
metric measurements of polar motion 


IVS station code for Richmond, FL 
root mean square 


Solution Independent Exchange 

satellite laser ranging 

a combination of independent space-geodetic measurements of the Earth’s 
orientation 

IVS station code for Svetloe, Russia 


transverse component of Earth orientation from single-baseline VLBI 
IVS station code for 32-m antenna at Tsukuba, Japan 


United States Naval Observatory 

Universal Time 

Universal Time as determined at a particular observatory 

The principle form of Universal Time; mean solar time at 0° longitude ob- 
tained by correcting UTO for observatory location and polar motion. 
Coordinated Universal Time 

University of Texas Center for Space Research 


vertical component of Earth orientation from single-baseline VLBI 
very long baseline interferometry 


variation of latitude 


IVS station code for Westford, MA 
IVS station code for Wettzell, Germany 
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